Magnetisation process in Er2Ti207 and Tb2Ti207 at 
very low temperature 



p. Bonvillefl§|, S. Petitj, I. Mirebeauf, J. Robertj, E. Lhotel§§ 
and C. Paulsen§§ 

f CEA, Centre de Saclay, Laboratoire Leon Brillouin, 91191 Gif-sur-Yvette, France 
t CEA, Centre de Saclay, IRAMIS/Service de Physique de I'Etat Condense, 91191 
Gif-sur-Yvette, France 

§§ Institut Neel, CNRS & Universite Joseph Fourier, 38042 Grenoble, France 

Abstract. We present a model which accounts for the high field magnetisation at 
very low temperature in two pyrochlore frustrated systems, Er2Ti207 and Tb2Ti207. 
The two compounds present very different ground states: Er2Ti207, which has a planar 
crystal field anisotropy, is an antiferromagnet with rN=l-2K, whereas Tb2Ti207 is 
expected to have Ising character and shows no magnetic ordering down to 0.05 K, 
being thus labelled a "spin liquid" . Our model is a mean field self-consistent calculation 
involving the 4 rare earth sites of a tetrahedron, the building unit of the pyrochlore 
lattice. It includes the full crystal field hamiltonian, the infinite range dipolar 
interaction and anisotropic nearest neighbour exchange described by a 3-component 
tensor. For Er2Ti207, we discuss the equivalence of our treatment of the exchange 
tensor, taken to be diagonal in a frame linked to a rare earth - rare earth bond, with 
the pseudo-spin hamiltonian recently developped for Kramers doublets in a pyrochlore 
lattice. In Tb2Ti207, an essential ingredient of our model is a symmetry breaking 
developping at very low temperature. We compare its prediction for the isothermal 
magnetisation with that of "the quantum spin ice" model. 
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1. Introduction 



The pyrochlore titanates, with formula R2Ti207 where R is a rare earth, have been the 
subject of intense studies for a decade [Ij. The pyrochlore lattice, made of tetrahedra 
joined by their vertices, leads indeed to a frustration of the exchange interaction in some 
specific situations. Owing to different crystal field properties, to the Kramers or non- 
Kramers character, to the relative importance of the dipolar interaction with respect 
to exchange, these compounds display a great variety of low temperature behaviours, 
the best known being the "spin-ice" ground state occurring in Ho2Ti207 and Dy2Ti207 
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E]. In the last years, there has been a break-through towards a full understanding of 
the low temperature properties of the Tb, Er and Yb members of the series, in particular 
with the recognition that the anisotropy of the nearest neighbour exchange interaction 
plays an essential role [U O |6]. We are interested here in the very low temperature field 
variation of the magnetisation in Er2Ti207 and Tb2Ti207, and we will show that much 
of the experimental data available in the two compounds can be accounted for, in each 
case, using a unique anisotropic exchange tensor. 

Er2Ti207 is an antiferromagnet with Tn=1.2K [7] and it presents an easy magnetic 
plane, perpendicular to the local < 111 > ternary axis, arising from the crystal field 
ground doublet of the Kramers Er^'^ ion |H [8]. The magnetic structure has been 
determined [9l [10] and magnetic order was suggested to arise from the "order by 
disorder" mechanism |Tl], which was recently put on a more robust ground [121 IT3] - 
The phase diagram and the evolution of the magnetic structure upon application of 
a magnetic field were also determined [HI [15] and interpreted in terms of anisotropic 
exchange, using a simple model [I5]. A symmetry constrained 4-component anisotropic 
exchange tensor was eventually introduced [16] and derived from fitting of the spin 
wave dispersion laws [I3]. In these latter works, the hamiltonian is written in terms of 
effective S=l/2 pseudo-spins which represent the ground doublet alone. This approach 
is adequate to describe the physics of Er2Ti207 at low temperature (the two first excited 
doublets have an energy of 73 K and 85 K above the ground state [9]) and in zero or 
moderate magnetic fields (a few T, of the order of the critical field He — 1.7T, see below), 
where mixing with the excited doublet is negligible. The pseudo-spin hamiltonian cannot 
however capture the high field magnetisation since mixing plays herein an important role. 
In the first part of this work, we start by showing the equivalence of the 4-component 
pseudo-spin exchange tensor and of the exchange tensor used in the present work, which 
is diagonal in a R-R bond frame. Then we show that the zero-field antiferromagnetic 
(AF) ground state in Er2Ti207 and the single crystal magnetisation curves in the AF 
phase are correctly described by mean field theory, considering the full crystal field 
interaction, involving thus the total rare earth momentum J, and anisotropic exchange. 
We show that the same exchange tensor also allows to reproduce the neutron diffuse 
scattering in the paramagnetic phase. 

In Tb2Ti207, a quite different picture is relevant. The ground state presents no 
long range magnetic ordering down to very low temperature [T7]. The crystal field level 
scheme of the non-Kramers Tb'^^ ion has two ground doublets separated by a small 
gap A ~1.4meV [13 [19], each having Ising anisotropy along the < 111 > ternary axis. 
The exchange in Tb2Ti207 is of AF type, and thus the lack of magnetic ordering is 
puzzling since antiferromagnetic exchange is not frustrated with Ising-like spins on the 
pyrochlore lattice. A model involving quantum fluctuations between the two ground 
doublets [20] was developed to account for the "spin-liquid" behaviour in Tb2Ti207. 
This model predicts that the very low temperature magnetisation curve should present 
a plateau (for a field along [111]) around 0.05 T [21], akin to the "2/3 plateau" observed 
in spin-ices [22]. However, two recent works [23t l2l] did not confirm this prediction. An 
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alternative model for Tb2Ti207 was proposed by us based on the occurence of a 

symmetry breaking at low temperature, caused for instance by a Jahn- Teller distortion 
of the rare earth site or by quadrupolar ordering. In a certain region in the space of 
exchange integrals, this approach yields a phase with no long range magnetic ordering. 
In the second part of this work, we show that the very low temperature magnetisation 
curves in a single crystal [231 l2l] can be reproduced, to a good approximation, in the 
frame of our model, using the anisotropic exchange tensor derived in Ref. |26]. 



2. The exchange tensor and the magnetisation in Er2Ti207 
2.1. The anisotropic exchange tensor 

Various conventions have been used to define the anisotropic exchange tensor in the 
pyrochlores. Since exchange is a two-ion interaction, a natural choice is to consider the 
vector linking two rare earth ions (the R-R bond) as principal axis for the exchange 
tensor |27]. We make this choice in the present work, i.e. we take a 3-component 
exchange tensor J' which is diagonal in a "bond frame" (a,b,c) where c is the R-R 
bond axis (see the Appendix for a definition of the bond frames). One can also add an 
antisymmetric Dzyaloshinsky-Moriya exchange term, with parameter J4. In this frame, 
the exchange interaction between ions i and j writes: 

where the are the components of the full angular momentum. Another choice is made 
in Refs. [T3l[T6] . where an S=l/2 pseudo-spin is considered to describe the ground doublet 
and the exchange hamiltonian is written in terms of the pseudo-spin components. In 
case of a Kramers doublet in the local frame with trigonal symmetry, the 3-component 
axial g-tensor writes g = {g±, g±, gz}, and the projection onto the ground doublet yields 
the relationship: gjj = g S, where gj is the ionic Lande factor. 

In the Appendix, we demonstrate that these two formulations of the exchange 
hamiltonian are equivalent and we show how to derive the relationships between 
the {J7a, J7b, i7c, J4} tensor used here and the exchange parameters {Jzz, J±, J±±, Jz±} 
entering the exchange hamiltonian in terms of pseudo-spins. They write as follows (for 
a Kramers doublet): 

r9z,2 -Ja + 2Jc + 4J4 



J: 



ZZ 



J /S'X N 2 2j7a — 3j76 — jT'c + 4 J4 , , 



^9] 3 
± -^2 2j7a 3^ 
J 

9z9± Ja^ Jc~ -h 

3^2 



The reciprocal relations write: 



Ja = Jzz + 4(^)^(J± - J±±) - 4V^^ J.±] 

3 gz 9± 9z9± 
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We emphasize that the physical parameters here are J'a, Jib, Jc-, Ja, and the two 
components and g± of the g-tensor, the pseudo-spin constants being effective 
parameters. 

2.2. Isothermal magnetisation in the AF phase and diffuse scattering in the 
paramagnetic phase 

The magnetisation vs field has been measured in the AF phase of an Er2Ti207 single 
crystal, at 0.15 K, using a SQUID magnetometer equipped with a dilution refrigerator 
developed at the Institut Neel-CNRS. The magnetic field was applied along the 3 
symmetry directions [111], [HO] and [100], and the data are shown as black dots in 
Fig{Tl The demagnetisation factor was negligible for the 3 directions. For H // [111] 
and [110] (FigslU a and b), one observes a steep increase of the magnetisation with 
a slight upwards curvature up to a critical field ^1.6T; for H // [100] (Fig|T]c), 
the initial increase is linear and the critical field slightly higher: He ~1.8T. Above 
He, a much slower linear increase is observed up to the maximum field of 7T. We 
have computed the field variation of the magnetic structure in each case, using a self- 
consistent mean field calculation which involves the 4 Er sites of a tetrahedron, each 
being exchange coupled to its nearest neighbours, and taking into account the infinite 
range dipole-dipole interaction. The underlying asssumption of this type of calculation 
is that the magnetic structure under field has a k=0 propagation vector. This is the case 
when the field is applied along the [110] direction [HI [15] and it is probably so for H // 
[111] and [100], as can be concluded a posteriori from our study. The Er^"*" ion (J=15/2, 
gj=6/5) is described by its full angular momentum J submitted to a trigonal crystal 
field interaction with parameters as in Ref.pj. The g-tensor of the ground doublet has 
components: gz=2.6 and 5'_l=6.8 

Starting from the pseudo-spin exchange parameters derived for Er2Ti207 in Ref. [13] , 
relations (jS]) allow us to obtain the exchange integrals J7a, J/j,, J^c and J4. Actually, the 
set of pseudo-spin parameters {Jzz, J±, J±±, Jz±} includes the dipolar interaction limited 
to first neighbours, whereas our calculation makes use of the exchange only integrals 
and of the infinite range dipolar interaction. Then one must replace J7a, J^b and J7c in 
relations ([3]) respectively by Ja — D, Jh — D and Jc + 2D (see Appendix), where D is the 
characteristic first neighbour dipolar energy worth 0.022 K in Er2Ti207. One obtains, 
taking into account the error bars given in Ref. [T3] for the pseudo-spin parameters: 
Ja= 0.070 ± 0.048 K, = -0.055 ±0.01 K, = -0.075 ±0.060 K and J4 = -0.018 
± 0.015 K. Except for jTj,, the range of acceptable values is quite large. The isothermal 
magnetisation being quite sensitive to the particular values of the exchange parameters. 
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Figure 1. Isothermal magnetisation curves in Er2Ti207 close to 0.15 K, with the 
magnetic field applied along [111] (a), [110] (b) and [100] (c): the black dots are the 
experimental data and the solid line the calculated curve. Insets: calculated moduli 
of the Er moments for the 4 sites in a tetrahedron. For H // [111], the black curve 
corresponds to the Er site with its ternary axis parallel to [111], the other colors to the 
3 other sites. For H // [110], the black and red curves correspond to the a-sites, for 

which the local axis is at an angle of acos(y^) ~ 35.3° from the field, and the blue 
curve to the /3-sites, for which the local axis is perpendicular to the field. For H // 
[100], the 4 moduli have the same field variation. All the simulations were performed 
with the anisotropic exchange tensor J7a — 0.030 K, Ji, = — 0.050 K, J7c — — 0.105 K 
and J4—O. 

the calculation of its field variation at 0.15 K for H // [111], [110] and [100] allows us 
to precise their values and to reduce significantly the error bars. We obtain the best 
match to the data (see FigJT]) with the following parameter values: J7a= 0.030 ±0.005 K, 
= -0.050 ±0.005 K, = -0.105 ±0.01 K and J4 = ±0.005 K. The slight upward 
curvature at low field for H // [111] and [110] and its absence for H // [100], the critical 
field values and the linear increase of the magnetisation above Hf. are well reproduced. 
The derived exchange parameter values, while lying within the range determined in 
Ref. [13] from the fit of the spin wave dispersion laws, represent a much more precise set 
for the exchange tensor. 

We note that, with these exchange integral values, the ground configuration in zero 
field is found to be the so-called antiferromagnetic 1^2 state, as determined from neutron 
diffraction [9l[l0j. This result, which is obtained using mean field theory alone including 
the total crystal field interaction, had been evoked in Ref. [28]. However, we find that 
the Er moments in the ^2 state do not lie exactly in their easy plane, but have a small 
(~ 0.01 /xb) out-of-plane z-component. The associated mean field Neel temperature is 
Tn= 2.62 K, higher than the experimental value 1.2 K. This enhancement of the mean 
field Tn value with respect to the actual value is attributed in Ref. [13] to the effect of 
classical spin fluctuations. 

The calculated field variations of the moduli of the 4 Er moments in a tetrahedron 
are represented in the insets of Fig HI For H // [111], the Er moment at the site with its 
ternary axis parallel to the field has a rather simple behaviour. Its modulus (black curve 
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Figure 2. In Er2Ti207 at 0.3 K, with the magnetic field apphed along [110]: a) 
moduli of the Er moments for the 4 sites; b) angles of the 4 Er moments with their 
local ternary axis <111>; c) angles of the 4 Er moments with the field H. The data 
are from Ref. |15| and the simulations (solid lines) were performed with the anisotropic 
exchange tensor Ja = 0.030 K, Jb = -0.050 K, Jc = -0.105 K and Ja ^ 0. The red 
and black data points and lines correspond to a-sites, the blue data points and lines 
to the two (3 sites, which have the same behaviour as a function of field. Note in c) 
that the moment at the a-site with local axis [111] (black line) reverses its direction 
near the critical field. 



in the inset of FigUJa) is seen to strongly decrease as H increases, signalling a departure 
from the easy plane; it reaches a minimum value near the critical field ^1.6 T, and 
increases linearly with H above He. For H // [110], a similar behaviour is obtained for 
the a-sites (black and red curve in the inset of FigiUb). For both these field directions, 
there are three different behaviours for the moments as a function of field, which implies 
that the field direction is not a symmetry axis for the magnetic structure. By contrast, 
for H // [100] (Fig{T]c), a single behaviour is obtained for the 4 sites, suggesting that 
[100] is a symmetry axis for the field induced structure, which could be checked by 
in-field neutron diffraction. 

Neutron diffraction was performed at 0.3 K for H // [110] in Er2Ti207 [lll|T5] and 
the data were interpreted in Ref. [15] using, as a crude approximation, a two-component 
exchange tensor diagonal in the local frame. We replot here the data of Ref. [T5] (Figj2]) 
and show that the evolution of the magnetic structure can be reproduced using the 
presently derived 3-component exchange tensor diagonal in the bond frames (neglecting 
the Dzyaloshinski-Moriya exchange term). This gives a physical insight about the field 
evolution of the magnetic structures, which present similar features whatever the field 
direction. On Figj2] c, which depicts the field variation of the angles between the Er 
moments and the field, it is clear that the moments rotate towards the field direction up 
to the critical field value, then align along (or close to) the field. The critical field can 
thus be considered as the "spin flip" field of the AF structure, and the linear increase of 
the moment moduli (and of the magnetisation) above is caused by quantum mixing 
with the excited crystal field states. The calculated values of the angles are in good 
agreement with the data (Fig|2]b and c), but the calculated moment moduli above 
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Figure 3. Diffuse neutron scattering at 2 K in Er2Ti207 in the {hkO) plane (upper 
panel) and in the (hkk) plane (lower panel) of the reciprocal space: Left: experimental 
data from Ref. j29| : Right: simulated maps (at 3.2 K) using the exchange tensor 
Ja =0.030 K, Jb = -0.050 K and Jc = -0.105 K in the 4-site RPA approximation. 



(Figj2]a) are somehow overestimated. 

The neutron diffuse scattering in the paramagnetic phase at 2K reported in 
Ref. [29] is displayed in the left part of FiglHl It shows intense spots near the Bragg 
positions (111), (133), (311), (333) and less intense maxima near (000), (022) ... 
Our calculations (Figj3] right part) of this diffuse scattering are performed using the 
RPA approximation [30] with the anisotropic exchange tensor derived above from the 
magnetisation curve. The temperature was chosen at 3.2 K since the mean field Neel 
temperature with these parameters is 2.62 K. One can see that our simulations reproduce 
satisfactorily the experimental data of Ref. |29j. We note that they capture the high 
intensity diffuse spots at Bragg positions, precursor to the long range magnetic ordering 
with k = occuring at lower temperature, whereas the simulations performed in Ref. [29] 
do not. We believe this is due to the fact that the approach used in this latter work, i.e. 
the diagonalisation of the 4-site exchange/ dipolar hamiltonian on a tetrahedron, cannot 
describe the critical correlations above responsible for these diffuse Bragg spots. 
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3. Very low temperature isothermal magnetisation in Tb2Ti207. 

In the spin liquid Tb2Ti207, where long range magnetic ordering does not occur down 
to very low temperature in zero field, application of a magnetic field induces magnetic 
order with k=0 for H // [110] |31] and H // [111] [32], and probably also for H // [001] 
although no neutron diffraction data are available for this field direction. Actually, for 
H // [110], an AF structure with k = [100] coexists with the k=0 structure above 2T 
and below 1 K |31], but it is not expected to contribute to the magnetisation. Then, the 
4-site self-consistent calculation sketched in the previous section is expected to hold for 
describing the magnetisation associated with the field induced magnetic structure, but 
a knowledge of the zero field ground state in Tb2Ti207 is necessary for this purpose. We 
recently proposed a model which can account for the zero-field spin liquid phase [251 126] , 
and which involves the development, at low temperature, of a symmetry breaking at the 
rare earth site. This symmetry breaking could be either an effect precursor to a Jahn- 
Teller transition occuring at much lower temperature or due to quadrupolar ordering, 
as recently suggested for non-Kramers ions [33]. The associated distortion, assumed to 
be of tetragonal symmetry [26l [33] , reads in the local frame, taking for instance its axis 
along the cubic [001] axis: 

•Hq = ^ [2J2 + J2 + (4J, + J,J,)], (4) 

where Dq is the strength of the distortion. The main effect of the symmetry breaking 
is to lift the degeneracy of the ground crystal field doublet and to destroy the Ising 
character of its wave-functions. It also results in the appearance, in the space of 
AF exchange integrals {J^a,<Jb,<Jc}, of regions where the mean field Tb moment is 
zero near zero temperature, i.e. of a spin liquid phase where short range correlations 
alone are present. The distortion strength and the exchange tensor derived in Ref . |26j . 
which describe correctly the temperature variation of the local susceptibility [1], the 
variation of the field-induced magnetic structure [31] and the inelastic and diffuse 
neutron scattering at very low temperature [31] in Tb2Ti207, are respectively: Dq 
= 0.25 K and Ja = -0.07K, Jb = -0.19 K and Jc = -0.09 K. In a single crystalline 
sample, one expects the tetragonal distortion to be distributed in domains with axis 
along the 3 fourfold cubic axes. Another feature in Tb2Ti207 is the magneto-elastic 
(ME) interaction which yields giant magnetostriction effects, due to the presence of 
the low lying crystal field level [351 EH]. Averaging over the domains and taking into 
account the magneto-elastic interaction (limited to quadratric terms in the total angular 
momentum) following the formalism of Ref. [36], we have computed the 0.08 K and 4K 
magnetisation for fields along the 3 symmetry directions [100], [110] and [111]. We used 
the above quoted anisotropic exchange tensor and the trigonal crystal field as in Ref. [4] , 
in the presence of the symmetry breaking (j4]). At 4 K, we set the distortion to zero since, 
at this temperature, its strength should be much smaller, but probably non- vanishing. 

The comparison between experimental data and our calculation is shown in FigjH 
The overall agreement is reasonably good for both temperatures 0.08 K and 4 K. In order 
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Figure 4. Isothermal magnetisation curves in Tb2Ti207 at 0.08 K and 4K: a for H // 
[111], b for H // [100] and [110]. The experimental data for H // [111] and [110] are 
taken from Ref. [U . The calculated curves were obtained using our model (tetragonal 
distortion and anisotropic exchange) with two assumptions: no magneto-elastic effects 
(dashed lines) and including the magneto-elastic distortions (solid lines), except for the 
solid blue line in a which reproduces the prediction of the "quantum spin ice" model 
at O.IK of Ref. ^21^, up to the highest field of the calculation (0.5 T). 

to assess the importance of the ME interactions, we have represented the calculated 
curves obtained without (dashed lines) and with (solid lines) ME effects. Inclusion of 
the ME interaction slightly modifies the magnetisation, especially at high fields. It 
yields a small enhancement, which results in a better agreement with experiment for 
H // [100] and [110] (FigJlb), but not for H // [111] (FigJla). For this latter field 
direction, we have also reproduced the curve calculated using the "quantum spin ice" 
model at 0.1 K (blue line), taken from Ref. j21j. Despite the limited field range, it is clear 
that it does not reproduce the data, which invalidates the "quantum spin ice" model for 
Tb2Ti207, at least as far as the magnetisation is concerned. Furthermore, this model 
predicts a sizeable variation of the shape of the low field magnetisation curve between 
0.02 K and O.IK, with the appearance of a clear plateau at the lowest temperatures 
(0.02 K). By contrast, the experimental data [23l [24] (and our calculations) show that 
the shape of the magnetisation curve does not appreciably change between 0.05 K and 
0.3 K. We believe the absence in Tb2Ti207 of the magnetisation plateau expected for 
Ising spins is due to the proposed symmetry breaking at the rare earth site, which 
destroys the Ising character linked with the bare trigonal crystal field wave-functions. 

Although our model correctly reproduces the overall magnetisation behaviour and 
its anisotropy, some deviation occurs at 0.08 K at low field, around 1 T and below, mainly 
for H // [111] and [100]. The curvature of the magnetisation as the field increases is not 
exactly reproduced. At higher fields, above 3-4 T, the calculated points lie somewhat 
above the data points, especially for H // [111] both at 0.08 K and 4 K. Keeping in mind 
that the uncertainty for such magnetic measurements is usually estimated to amount 
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to a few percent, we can envisage various causes for these deviations. First, our model 
would not capture all the details of the field-induced magnetic structure, especially at 
low fields. Second, the effect of a slight field misalignment with respect to the crystal 
axes can also play a role. For H // [110], indeed, there occurs a "spin melting" near 1 T 
where the two moments lying on sites with their ternary axis perpendicular to the field 
[P sites) vanish [31]. The occurence of the "spin melting" and the configuration of the P 
moments are very sensitive to the alignment of the field with respect to the crystal axis. 
For H // [111], a "spin melting" seems also to occur, but only for one of the 4 Tb sites and 
at a much lower field [32], and a small misalignment of the field should have practically 
no effect on the computed curves. Finally, we made an approximation concerning the 
tetragonal distortion and the magneto-elastic interaction. We assumed that they can 
be treated independently, i.e. that the field induced strains merely superimpose onto 
the tetragonal field independent strain. This implies that the distortion domains are 
equiprobable and that they are not affected by the magnetic field. At low field, where 
the ME term in the hamiltonian has a magnitude of a few 10^^ K, to be compared with 
Dq = 0.25 K, this assumption seems valid, but at high field where the ME interaction 
can reach a few O.IK, it could be questionable. A more elaborate theory would be 
needed to account for these effects, including eventually higher order terms and/or 
other interactions like the quadrupole-quadrupole coupling, and their field dependence. 

4. Conclusion 

We have presented a mean field approach for frustrated pyrochlore systems which takes 
into account the full crystal field level scheme of the rare earth ion, the anisotropic 
nearest neighbour exchange and the infinite range dipolar interaction. The model 
is applied to investigate the high field low temperature magnetisation curves in the 
antiferromagnet Er2Ti207 and the spin liquid Tb2Ti207, together with other physical 
properties. In Er2Ti207, we show the equivalence of the exchange tensor diagonal in the 
bond frame we use here and of the pseudo-spin S=l/2 exchange tensor recently proposed. 
Our fit of the magnetisation curve in the AF phase, where the quantum mixing with 
excited crystal field states is well accounted for, allows us to derive much more precise 
values for the exchange parameters. These parameters also describe well the diffuse 
neutron scattering in the paramagnetic phase. In Tb2Ti207, we introduced a modified 
crystal field with a small tetragonal distortion from trigonal symmetry and we took into 
account the magneto-elastic interaction. We show that the exchange tensor we derived in 
previous works allows the low temperature magnetisation curves for high symmetry field 
directions to be reasonably well reproduced without further parameters. The advantage 
of the present approach with respect to models dealing with S=l/2 pseudo-spins lies in 
the fact that it can take into account the effects linked with quantum mixing of crystal 
field states by the magnetic field. This is important in Er2Ti207 when dealing with 
high field properties, and still more important in Tb2Ti207 since the first crystal field 
excitation has an energy of ~15K and mixing is important even for moderate magnetic 
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fields. 
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Appendix A. Relations between the exchange tensor appropriate to the 
S=l/2 pseudo-spin hamiltonian and the exchange tensor in the "bond 
frame" 

We start by choosing a set of vectors forming an orthonormal frame attached to the 
bond linking rare earth neighbours i and j. Calling the unit vector along the < 111 > 
trigonal axis at site i, we define the following unit vectors forming the "bond frame", 
where lies along the link between the two neighboring sites : 

au=^(ei + ej) (A.l) 
b^j c^j X a^j . 

In the bond frame, we consider the general anisotropic exchange Hamiltonian, written 
in terms of the full angular momenta J and where the last term is the antisymmetric 
Dzyaloshinski-Moriya exchange (the convention used here is that an antiferromagnetic 
exchange integral is negative) : 

^ = -[ X] Ja(Ji-aij-)(aij-Jj) + J6(Ji-bij)(bij-.Jj) 

<ij> 

+ j;(Ji.Cij-)(cii-Jj) + J4 bij.Ji X Jj], (A.2) 

where Z]<jj> means a summation over the first neighbour pairs. This can be written in 
terms of a global exchange matrix J: 

^ = - E Jii'^a al al + J, hi hi + Je cl cl + V2 J4 nuv) J] 

<ij>,uv 

= - ^ J,J3„ (A.3) 

<ij> 

where fj is an antisymmetric 3x3 matrix with only two non-zero elements: 7713=1 and 
?73i = — 1. In the "bond frame", the matrix S writes: 

/ Ja V2JA 



J 



Jb 
{ -V2 J4 J 



(A.4) 
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In Ref . [13] , the general anisotropic exchange Hamiltonian is written in terms of the 
spin components within the local frame of an S=l/2 pseudo-spin: 



(A.5) 



where the "sanserif notations refer to spin components in the local bases. The two 
states of this S=l/2 pseudo-spin span the states of the ground CEF doublet and, for 
the case of a Kramers ion (like Er^"^), one can define a local g-tensor g such that: 
gjj = gS. In the context of pyrochlores with local trigonal symmetry, the g matrix is 
diagonal and takes the form: 



9 



( 



9± 







9± 




\ 


9z J 



(A.6) 



Our goal is to determine the relation between the "bond frame" exchange 
parameters J'a, Jh-, Jc and J4, on the one hand, and the pseudo-spin parameters J^^, 
Jzit, J± and J±± on the other hand. To this end, we call Bij (resp. Mi) the matrix 
transforming the coordinates in the "bond frame" {ij} (resp. the local frame i) to the 
cubic (cartesian) coordinates, and A the matrix transforming (Sx, Sy, Sz) into (S+, S_, Sz) 
in the local basis (we omitted the site indexes for sake of clarity): (?jJ = g A with 
S = (S+, S_, Sz), and: 



A 



V 



1/2 l/2i \ 
-l/2z 1/2 

1 y 



(A.7) 



A straightforward transformation of hamiltonian ( 1A.3P shows that: 



n 




T ~ 



g Mj B,^ J Bf^ ~g A 



(A.8) 



This allows the relations between the two sets of parameters appearing in ( 1A.4I) and ( lA.Sp 
to be determined; they are given in section |2TT] of the main text. The conventions as to 
the sign of the exchange integrals are different for the two sets: for an AF interaction, 
they are positive (except J±) in the pseudo-spin representation and negative in the 
"bond frame" description. 

We note that, in the "bond frame", the dipolar interaction limited to a nearest 
neighbor pair {ij}: H^iip = D [JjJj — 3 (Jj.Cjj) (cjj.Jj)] = 3i D 3j is diagonal with: 



and D 



D 



4-71 fce 



''dip 

( D 


V 




D 







-2D 



\ 



(A.9) 



^^5^, where a is the cubic cell lattice parameter. 
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